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We develop a theory of continuous decoupling with bounded controls from a geometric perspective. 
Continuous decoupling with bounded controls can accomplish the same decoupling effect as the bang- 
bang control while using realistic control resources and it is robust against systematic implementation 
errors. We show that the decoupling condition within this framework is equivalent to average out 
error vectors whose trajectories are determined by the control Hamiltonian. The decoupling pulses 
can be intuitively designed once the structure function of the corresponding SU (n) is known and is 
represented from the geometric perspective. Several examples are given to illustrate the basic idea. 
From the physical implementation point of view we argue that the efficiency of the decoupling is 
determined not by the order of the decoupling group but by the minimal time required to finish a 
decoupling cycle. 



I. INTRODUCTION 

Quantum decoherence control has been studied inten- 
sively in recent years due to the fact that it represents 
one of the main obstacles in implementing quantum com- 
putation. However the problem of decoherence remains 
daunting. Although there have been rapid advances in 
physical realization of quantum operations on few lev- 
els system in which a single qubit or few qubit can be 
embedded . The fidelity for those quantum operations, 
such as Rabi rotation, is typically far below the thresh- 
old to make the system scalable. To increase the fidelity, 
many strategies are proposed to counteract the unde- 
sired effects resulting from decoherence. Those strategies 
can be roughly divided into two categories, depended on 
if an encoding scheme is used. Some prominent exam- 
ples in which encoding is needed are the quantum error- 
correction codes (QECCQElllOl and the decoherence- 
free subspace (DFS).0,Hlg One of the main drawback 
of the encoding strategies is the large amount of ancil- 
lary space needed, especially when the qubit is still a rare 
resource. As a result, strategies which do not need an- 
cillary systems is highly desirable. Closed-loop decoher- 
ence control (such as quantum feedback control 
and open-loop decoherence control (such as bang-bang 
or dynamical decoupling |ToL ITTl IT^ I fall into this 
category. The open-loop decoherence control typically 
requires only a limited, a priori, knowledge of the system- 
environment dynamics. Empirical determination of con- 
trol parameters via quantum tomography has also been 
proposed recently. |l4| In this work we will focus on the 
dynamical decoupling, but it is acknowledged that no 
single strategy can efficiently eliminated the problem of 
decoherence at all levels. Combining or concatenating 
different strategies is usually necessary for any real phys- 
ical system. 

In the original bang-bang decoupling framework, arbi- 
trarily strong and instantaneous control pulses are uti- 
lized to induce frequent unitary interruptions during the 
evolution of the system. The control Hamiltonian is in- 
dependent of the system dynamics and are judicially de- 
signed to realize an effective decoupling between the sys- 



tem and the environment. The stringent requirement 
on the control Hamiltonian represents a drawback of the 
bang-bang decoupling. Experimentally it is impossible 
to implement arbitrarily strong and instantaneous pulses 
in real physical systems. An ultra strong and fast control 
pulse will also inevitably induce transitions to higher en- 
ergy levels which are usually neglected in the analysis of 
the bang-bang decoupling. Theoretically is difficulty to 
describe the evolution with and without control terms si- 
multaneously. This also make it difficult to estimate the 
robustness of the control and to estimate the error in- 
duced by operational imperfections. The highly abstract 
group theory which is frequently used in the decoupling 
analysis also make it less intuitive and difficult to make 
transparent connection to the realization of bang-bang 
decoupling. To alleviate the need for strong and impul- 
sive control actions, a dynamical decoupling using only 
bounded-strength Hamiltonian is recently proposed, [lq 
Within this framework the same group symmetrization 
is achieved by exploiting the Eulerian cycles on a Cay- 
ley graph of £7, [16] where Q is the decoupling group. On 
the other hand a complimentary, geometric perspective 
of bang-bang decoupling is recently proposed to provide 
a more intuitive picture and provide a method to esti- 
mate implementation errors. |l4l Il7| In this work we ex- 
tend the geometric picture of bang-bang decoupling to 
the case of continuous decoupling with bounded controls. 
We show that the decoupling by symmetrization is equiv- 
alent to average out error vectors whose trajectories are 
determined by the control Hamiltonians. The decoupling 
pulses can be intuitively designed by viewing the effect 
of control Hamiltonian from a geometric perspective. To 
avoid confusing we will use the term bang-bang decou- 
pling to refer to the typical decoupling scheme originally 
proposed. |ld ITU Il2t Il3 | The term Eulerian decoupling 
corresponds to the decoupling with bounded Hamiltonian 
introduced by L. Viola and E. Knill^f|, while the term 
continuous decoupling corresponds to the version of con- 
tinuous decoupling from a geometric picture developed 
in this work. 

The structure of the paper is the following. In Sec. 
[n]we review the original bang-bang decoupling, the Eu- 
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lerian decoupling, and the geometric perspective of the 
bang-bang decoupling. In Sec. 11111 we introduce the 
idea of continuous dynamical decoupling from a geomet- 
ric perspective. In Sec. IIVI several examples are given 
to illustrate how to construct the decoupling pulses. We 
summarize in Sec. 



II. REVIEW 

A. Bang-bang decoupling decoupling 

In general, dynamical decoupling seeks to eliminate 
the decoherence of an open quantum system by effec- 
tively averaging out the interaction between the system 
and the bath by introducing some stro ng p eriodic control 
Hamiltonian on the system. 0, EH El El Following the 
standard treatment, the dynamics of an quantum open 
system is determined by the total Hamiltonian 



where 



H a =H S ® 1b + 1s®H b + H S b, (1) 



H SB = Y,S 1 ®B 1 (2) 



represents the system-bath interaction. Now introduce a 
control Hamiltonian H c (t) = H c (t) ® lg acting on the 
system alone. Denote U c (t) the time-evolution operator 
associated with the control Hamiltonian 



U c (t) 



Texp |- 



duH c (u) 



(3) 



where T is the time-ordering operator. In periodic dy- 
namical decoupling one is restricted to the situation 
where the control field is cyclic, i.e., U c (t) — U c {t + T c ) 
for some period T c . For any state \ip(t)) and any op- 
erator O in Schrodinger picture, the corresponding state 
\ip(t)) and operator 0(t) in the interaction representation 
associated with H c (t) are: 



\m) = umm) = utwwMo)), (4) 



and 



d{t) = u*(t)ou c (t), 



(5) 



where U(t) = Texp {—iHt} is the time-evolution oper- 
ator associated with H = H + H c . It is easy to show 
that 



\ip(t)} = Texp 



n: 



duH o ( U )\\<p(0)), (6) 



where Ho(t) — U^.(t)HoU c (t). Using the periodic condi- 
tion U C (T C ) — Is one finds 

\1>(T e )) = U(T c )\m) = Ut(T c )U(T c )\i>(0)) = \i,(T c )). 

(7) 



The stroboscopic dynamics at T/v = NT C for integer N is 
hence determined by the time evolution operator 



U(T C ) = Texp ■ 



-i I duHo(u) 



(8) 



It is possible to define a k-th order average Hamilto- 
nian H k by applying the Magnus expansion (18| to U(T C ), 
results in 



U(T C ) : 

In particular we have 



,-i[H°+B 1 +...]Tc 



(9) 



H v = 



1 



duHo(u) 



1 



duUi{u)H (u)U c {u), 



and 



H x = -, 



T 



dv I du 
lo 



H (v),H (u) 



(10) 



(11) 



Higher-order corrections can also be systematically eval- 
uated. In the limit of T c — * 0, which is the ideal limit of 
bang-bang decoupling, H° becomes the dominant term. 
It can be viewed as an effective Hamiltonian governing 
the stroboscopic dynamics under control Hamiltonian. 
In the following we will focus on designing the effective 
Hamiltonian H°. 

Denote Tts the finite-dimensional state space associ- 
ated with the system. The idea of decoupling by sym- 
metrization is to identify a discrete decoupling group 
S = {9j}, J = {0, ■ ■ ■ , \G\ - 1}, acting on H s via a 
faithful, unitary, projective representation fi such that 
fJ'iG) C U(TLs)i the unitary matrices acting on TLs- The 
bang-bang decoupling via Q is implemented by assigning 



U C {{1 - l)Af + s) = fx(gi), s g [0, At), 



(12) 



with T c = |^| At, and I = {0, . . . , \Q\- 1}. With this 
assignment the effective Hamiltonian H° becomes: 



where 



Sry — Tig [Sry 



E^ 1 

7 



B~, 



11 s 3 ee 



(13) 



(14) 



The nontrivial work is to identify the group Q such that 
for all 7 the effective error operator S 7 = A 7 ls where 
A 7 is a real number. Once this is accomplished then the 
effective Hamiltonian H° is reduced to 



H° = lc 



(15) 



As a result the system is effectively decoupled from the 
bath, or equivalently the decoherence is suppressed. Note 
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that in this formulation the underlying control Hamilto- 
nian is never explicitly mentioned. An instantaneous, 
arbitrary strong control Hamiltonian is needed to imple- 
ment the desired U c (t). However, physically it is impos- 
sible to implement such an instantaneous control pulse. 
An arbitrary strong control Hamiltonian would also in- 
evitably induce transition to higher energy states which 
are neglected when a two level approximation is used 
to describe the qubit space. Those unphysical require- 
ments represent some of the main drawbacks of the orig- 
inal bang-bang decoupling framework. 



B. Eulerian decoupling 

In order to alleviate the unphysical requirements of 
the original bang-bang decoupling framework, L. Viola 
and E. Knill propose a general framework in which the 
same group symmetrization can be achieved while us- 
ing only bounded control Hamiltonians.^BJ Physically it 
uses only bounded control Hamiltonians to steer the time 
evolution operator. Mathematically it exploits the Eule- 
rian cycles on a Cayley graph^fj of the decoupling group 
Q . Given a decoupling group Q, the first step in imple- 
menting Eulerian decoupling is to find a generating set 
T = {f a }, a — 1) • • • ; 1-^*1) for the decoupling group Q. 
The physical implementation requirement is the ability 
to generate f a by some control Hamiltonian h a (t) over a 
period of time At, 



f a = Texp ■ 



At 1 

i I duh a (u) > , a = 1, . . . , [J 7 ]. (16) 



If we image each group element gi € Q as a vertex, 
then f a can be imaged as the directional, colored edge 
connecting the vertices. If g' = f a g, then we draw a line 
from point g to pint g' with color a. An Eulerian cycle 
is defined as a cycle that uses each edge exactly once. 
In this case, one can show that it is always possible to 
find Eulerian cycle, having length L = \G\\J~\. 16] A 
well-defined Eulerian cycle beginning at the identity go 
of Q can be uniquely specified by the sequence of the 
edge colors used, Pe = (pi,P2, ■ ■ ■ ,Pl), where pi € T . 
An Eulerian decoupling is then implemented by letting 
T c = LAt and by assigning U c {t) as follows: 



U c {t)[{l - l)At + s} = ui(s)U c [(l - l)At], 



(17) 



where s 6 [0, At), and ui(s) = Texp{~i J Q S duhi(u)}, 
ui(At) = fi(pi). In this way the average Hamiltonian 
H° becomes: 



(18) 



where 



\H5\AtJ 



l^l|ff|At 



duUt{t)H {t)U c {t) 



gj€G 



gmAt 



At 



dtu a (t)S 7 u a (t) > fi(gj 



It can be shown that the same decoupling can be 
achieved through this average. Assuming that At re- 
mains the same as in the Eulerian decoupling, the length 
of the Eulerian decoupling is lengthened by a factor of 
\F\ compared to the bang-bang decoupling. Eulerian de- 
coupling provides a guideline to design the control Hamil- 
tonian if the decoupling group and it's representation is 
known. However the search for the decoupling group re- 
mains a nontrivial work. 



C. The geometric perspective of the bang-bang 
decoupling 

Bang-bang decoupling and Eulerian decoupling 
schemes make heavy use of the abstract group theory. 
In recent years a complementary geometric perspective 
is developed to provide a more intuitive picture. ^3 
The geometric picture of the bang-bang decoupling 
utilities the Homomorphic mapping between the Lie 
group SU{n) and SO(N), where N = n 2 - 1. Let A», 
i = 1, . . . , N be the N traceless, Hermitian generators of 
SU(n). The generators satisfy trace-orthogonality, 



Tr(\i\j 



(19) 



where M is a normalization constant. For any group el- 
ement U € SU(n) 7 one can define a rotation R 6 SO(N) 
via 



A? 



(20) 



This defines a homomorphism from SU(n) to a subgroup 
of SO(N). 

Without loss of generality, one can re-write H as 

#o (21) 
= H s ® 1b + Is ® H B + H SB 

= Tv(H s )ts ® 1b + {H s - Tr(if<j)l) ® 1 B + l s ® H B 
+ l s ®^Tr(S' 7 )S 7 + ^(S' 7 -Tr( > S' 7 ))®S 7 

7 7 

= # + 1 S <8> #B + XX ^ B 7> 

7' 

where S' is traceless. Note that the first term Eq = 
Tr(i?g) only gives rise to an overall phase and can be 
discarded. For simplicity (and without loss of generality) 
we will write Hsb = E 7 Sy®B^, where S 1 is traceless. 

Any traceless system operator Sj € SU (n) can be ex- 
panded in terms of , yielding: 



S>y — ^ ^ (^7)2^2 — ^7 ' Xj 



(22) 
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where 



(* 7 )< = ]fJ Tr (A^ 7 ). 



(23) 



In other word, a traceless system operator can be repre- 
sented by a iV-dimensional vector . Using this result the 
system-bath Hamiltonian Hsb can be written as follows: 



H SB = E 5 7 ® B i = E & ■^)^ B i 



(24) 



where s 7 is a vector of length N. In the following we 
should refer s 7 as error vectors. By using error vectors 
to represent the system-bath Hamiltonian H$b, the de- 
coupling by symmetrization over a group Q with pro- 
jective representation fi(gi) can be viewed as an average 
over rotated error vectors. Following this line the average 
Hamiltonian H° becomes 



n ° = EiHiEM^Mff)®^ (25) 
7 ' y ' see 

= E i^i E^^EW'^)® B i 

7 11 see i 

= E T^j E E(*7)i-R[M<7)]y Aj ® S 7 
7 see y 

= EE 1 EE ^im^m^ [ ^ ® s 7 

7 j [ 11 see » J 

= (P 7 • X) <g> B 7 , 



where the average error vector s' 7 is equal to 



see 



see 



(26) 

From the geometric perspective the decoupling condition 
(in case of quantum memory) is equal to require that the 
average error s 7 be zero for all 7. Geometrically each 
term in group symmetrization procedure corresponds to 
an effective rotation R^[/i(g)] e O(N) on all error vec- 
tors. However it is evident that the error vector can be 
averaged to zero by a set of rotations which do not corre- 
spond to the representation of some decoupling group. It 
is thus intriguing to discuss if a underlying group struc- 
ture is necessary to achieve dynamical decoupling. 



III. CONTINUOUS DECOUPLING FROM A 
GEOMETRIC PERSPECTIVE 

The geometric picture of the bang-bang decoupling is 
intuitive but it shares the same drawback as the bang- 
bang decoupling, i.e., the error vector is instantaneously 
rotated to another vector by some rotation. This draw- 
back can be alleviated if a bounded control Hamiltonian 



is used to continuously rotate the error vector. We thus 
seek to formulate a framework for continuous dynamical 
decoupling from a geometric perspective. Recall that the 
average Hamiltonian has the following expression 



1 



l c Jo 



duU^(u)H (u)U c (u) =J2 S -y® B -y ( 27 ) 



Now instead of piecewisely mapping U c (t) into the rep- 
resentation of some decoupling group, one represents the 
effective system operator 5 7 by its corresponding aver- 
age error vector s' 7 . The average error vector can be 
expressed as the time average over the trajectory of the 
error vector rotated by [U c ] : 



1 f Tc 

L c Jo 



(28) 



1 / ' duY,(sM(u)XiU c (u) 
1 f Tc 



s' 7 • A, 



where 



L c Jo 



«); = 7f / C duJ2RHU c (u)Us,) t . (29) 



Or using vector notation: 



'-c Jo 



duR*[U c (u)]s^ 



(30) 



The decoupling condition (in case of quantum mem- 
ory) is to require s' 7 = for all 7. Note that we don't 
explicitly require that U c be the representation of some 
group. In order to make the design of control Hamil- 
tonian easier it is desirable to make a more transpar- 
ent connection between the decoupling and the control 
Hamiltonian. First define a time-dependent error vector 
via 

s- Ul(t)XU c (t) = s(t) ■ A. (31) 
The i-th component of s(t) can be expressed as 

Si(t) = ^Tr [\ t s- Ut{t)\U c (t)) . (32) 

It is instructive to study the trajectory of s(t) when the 
control Hamiltonian is proportional to one of the gen- 
erators of SU(n). Assuming that Hi(t) = ai(t)Xi, where 



■5 



ai(t) represents the envelope function of the control pulse, 
one finds 



d_ 

df 



M Si {t) = ^[Xiiai^hs-Ul^XUS)} 

+ Tr {\is ■ Ul{t)\U c {t){-m{t))h)} 



(33) 



i(n(t)Tr{ Xi 



Xi^sjtyXj 



= iai(t)y^Sj(t)Tr ^Xjiy^fijkXk j 



where we have defined 

[ L i)ij = -ifiji = +iflij, 



(34) 



and fijk is the structure function of SU(n). From the 
O(N) point of view, the effect of control Hamiltonian 
ai(t)Xi is to induce a rotation with generator Li and with 
speed ai (t) . It is thus useful to express Li in terms of the 
natural generators of SO(N). The natural generator of 
SO(N) are antisymmetric Hermitian matrices L^ u where 
[i, v = 1, • • • N, whose components have the form 



(35) 



We will restrict ourself to ^ < v as a convention and 
to avoid double counting. They satisfy the commutation 
relation 



-i - 5 



(36) 



/J,ft' L VV ' 



and the trace orthogonality condition 
Tr^L^,) = M5^5„ 



(37) 



Giving a control Hamiltonian of the form ai {t)Xi , in the 
geometric picture it corresponds to a time independent 
rotation generator Li and an time-dependent envelope 
function ai(t) representing the time dependent speed of 
the rotation. By expressing Li in terms of the natural 
generators of SO(N) 



a compact notation of the form 



(38) 



(39) 



can be used to represent the effect of control Hamiltonian 
A;. For example, for SU(3) — > SO(8) one has 



[A 3 ,A] = (-L12 + L 45 + L 76 ) X. 



(40) 



It immediately leads us to the conclusion that A3 can be 
used to average any vector in the 1-2 plan to zero with 
appropriate envelope function. In the appendix we ex- 
plicitly calculate the [A;, A] for SU(2), SU(3), and SU(4). 
Using the continuous decoupling from a geometric per- 
spective, the design of the decoupling pulse sequences 
can be outlined as follows: The first step is to find the 
corresponding error vectors from the system-bath Hamil- 
tonian. The second step is to identify the useful control 
Hamiltonian A;. The third step is to design a proper en- 
velope function a(t) to ensure the average error vector to 
be zero. 

Continuous decoupling is robust against implementa- 
tion imperfection because a small implementation er- 
ror on the control Hamiltonian will only result in a 
small deviation of the average error vectors from their 
ideal values. The implementation error can be evalu- 
ated via the distance between the ideal average error 
vector and the real average error vectors. Let s'a,^ be 
ideal average error vectors resulted from a perfect con- 
trol Hamiltonian while s' \ be the real average error vec- 
tors resulted from a imperfect control Hamiltonian. The 
Euclidean distance between two vectors d(s' \ ic i, s' \) = 



{s'\,id — s'x) ■ {s'xjd — s'\) can be used to quantify the 
implementation error. 



IV. EXAMPLES 

In this section we study several examples which illus- 
trate the basic idea of continuous decoupling from geo- 
metric perspective. 



A. Single qubit spin-flip decoherence 

Consider a single qubit with a single error operator, 
{S a } — {<7 Z }- In geometric picture it corresponds to a 
single error vector s = (0,0,1). Intuitively, performing 
any 2ir rotation in a plan containing the error vector s 
should average the error vector to zero. Using the results 
in the appendix it is easy to verify that one can choose 
the control Hamiltonian to be proportional to Ai(= a x ) 
or A2(= <Jy). From the geometric perspective, this corre- 
sponds to rotate the error vector in x-z or y-z plan using 
50(3) generator L13 or £23- If we were to choose a y and 
assume that H c (t) = a(t)<j y , the time-dependent error 
vector becomes 

s(t) = e -^l dua{u)Ll3 s(Q) = e- tA ^ L ^s(0), (41) 

where A(t) = J^dua(u). The decoupling condition can 
be written as 



L c jo 



— / due- lA(t)Ll3 s(0) =0. 



(42) 



The decoupling condition can be satisfied very generally 
by requiring that 

e -iA(T c /2)L 13 =i?[ 7r ^] ; (43) 

and 

A{T c /2 + t)=n + A(t),t = [0, T c /2), (44) 

where R[n,y] represents a 7r rotation around y-axis. 

Geometrically this corresponds to rotate the error vec- 
tor from s to — s at some speed controlled by a(t), and 
rotate it back to s with the same speed profile. The first 
condition ensures that the error vector is steered to — s 
at half time t — T c /2. The second condition ensures that 
the contribution from the second half cancels exactly the 
contribution from the first half, resulting in zero average 
error vector. 

It is known that the minimal decoupling group of bang- 
bang decoupling in this case is the group Z2 = {e,g} 
where g 2 = e. It is also referred as parity kick in the lit- 
erature. The corresponding representation can be chosen 
to be either {f/,(e) = l,(J,(g) = cr x } or {/x(e) = %n{g) = 
a y }. From the geometric picture the effect of /z(<?) is 
to rotate error vector s instantaneously to —a, ensuring 
the average error vector be zero. It is easy to show that 
a larger group C4 = {e, g 2 , g 3 , g 3 }, where g 4 = e, can 
achieve the same decoupling. The corresponding repre- 
sentation is {//(e) = l,n(g) = r,/i(g 2 ) = r 2 ,n(g 2 ) = r 3 } 
where r is the tt/2 rotation along x-direction. From the 
geometric picture this corresponds to rotate the error vec- 
tor by tt/2 at each kick. From the typical bang-bang 
decoupling point of view larger decoupling group repre- 
sents a less optimal decoupling scheme since more kicks 
are needed. On the other hand when viewed as the lim- 
iting case of the continuous decoupling, Z2 bang-bang 
decoupling corresponds to require a large amplitude a(t) 
for small t and turn a(t) off once the error vector s is 
steered to — s. A different a(t) can be similarly designed 
to reproduce the C4 bang-bang decoupling. However in 
any real physical implementation there is a upper limit 
for the strength of a(t). Hence there is a minimal time 
needed to finish one continuous decoupling cycle and the 
ideal limit of bang-bang decoupling is never reached. It is 
this minimal T c , not the order of the decoupling group, 
when compared to the decoherence time, indicates the 
efficiency of the decoupling scheme. 

B. Single qubit full decoherence 

Consider next a single qubit with all possible error 
operators, {S a } = {<J x ,<J v ,a z }. In geometric picture 
they correspond to three error vectors s\ = (1,0,0), 
,?2 = (0,1,0), and S3 = (0,0,1). Intuitively a 2tt rota- 
tion in a plan containing both the vectors s\ and S2 can 
average both error vectors to zero, but S3 will remain 
unchanged. However it is possible to design a sequence 
which average all three error vectors to zero. To see how 




Time (Arbitrary unit) 



FIG. 1: Trajectories of error vectors under the decoupling 
sequence described in example B. It is assumed that the en- 
velope function in each sub-period is the same. 



such a pulse sequence can be constructed, first recall that 
in previous example the vector Si(t) must reach — s*j dur- 
ing the decoupling operation. We thus seek to rotate s$ to 
— Si using alternating generators, in hope that the resid- 
ual errors will cancel each other when a full decoupling 
cycle is finished. Using this idea it is straightforward to 
construct and verify that a sequence of it rotations using 
following generators {Ai, A3, Ai, A3, A3, Ai, A3, Ai} can av- 
erage all three error vectors to zero. In FigQ]we plot the 
trajectories of three error vectors during the decoupling 
sequence. It is evident from the figure that a judicial en- 
velope function will ensure that all time integrals of error 
vectors are zero. 

The decoupling pulse sequence designed here is equiv- 
alent to the Eulerian decoupling prescribed in Ref IT5L 
However as pointed out in previous example the effi- 
ciency consideration should be based on the continuous 
decoupling framework. In Refll7l an example is given to 
demonstrate that a decoupling group may not be neces- 
sary to achieve the decoupling condition. From the geo- 
metric picture of the bang-bang decoupling it suffices to 
rotate instantaneously the error vectors to the vortices of 
a tetrahedron. From continuous decoupling point of view 
it is evident that for this kind of decoupling sequences the 
residual errors will accumulate when the error vectors are 
rotated from one vertex to another. It is difficult to can- 
cel these residuals systematically using a design similar 
to what has been done in this example. Even if it is ac- 
complished, it would represent a less optimal solution. 
We thus argue that a decoupling group is not necessary 
to achieve continuous dynamical decoupling but a de- 
coupling group usually provides a guideline to design the 
optimal continuous decoupling sequence. 
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C. Two qubits independent decoherence 

Next example consists of a two qubit system with in- 
dependent dephasing error operators. The system-bath 
Hamiltonian has the form: 

H SB =9i{a z ® l))®B 1 +g 2 {l®a z )®B 2 . (45) 

The Hamiltonian corresponds to two error vectors s\ = 
and s 2 — 52^6- Any rotation in SO(15) which can 
rotate both error vectors by 2ir should effectively decou- 
ple the system from the dephasing error. In many solid- 
state qubit, Heisenberg exchange is utilized to implement 
quantum operations. It is thus advantageous to design 
a decoupling sequence which is compatible with Heisen- 
berg exchange interaction. In terms of Pauli matrices the 
Heisenberg interaction is written as 

H ex — J&x ® <J X + J&y ® <J y + J(T Z (g) O z . (46) 

In other words the Heisenberg interaction corresponds 
to three vectors A7, An, and A15 from the geometric 
perspective. A decoupling sequence compatible with 
Heisenberg interaction should average error vectors to 
zero while leave those three vectors in tact. Using the re- 
sults in the appendix it is apparent that one should avoid 
using the generators Ai • • ■ Ae, Ag • • • A9, and A12 • • • A14 
which either can't average error vectors to zero or have 
undesirable side effects on Heisenberg interaction. This 
leaves us to the only choice of A7, An, and A15. It is then 
easy to verify that any one of those three generators can 
be used as the control Hamiltonian to average out the 
dephasing error with an appropriately designed envelope 
function a(t) while leaves the Heisenberg interaction in 
tact. 



D. One qubit with one ancilla level 

Our final example consists of a three level system in 
which the first two levels |1), |2) are used as the qubit 
space while the third level |3) is used as an ancilla level 
used to implement quantum operation in qubit space. 
We assume that all three levels are coupled to the envi- 
ronment, with the following system-bath Hamiltonian: 

H SB = 5 i(|l)(l|-|2)(2|)®5 1 (47) 
+ g 2 (|1)(3| + |3)(I|) ® B 2 + .93 (|2)(3| + |3)(2|) ® B 3 . 

This system-bath Hamiltonian corresponds to three error 
vectors si = A3, s 2 = A4, and S3 = Xq. By consulting the 
SU(2) — > 50(8) mapping in the appendix it is easy to 
show that all three error vectors can be averaged to zero 
via the generator A2 . This simple example illustrates that 
the continuous dynamical decoupling from geometric per- 
spective can be straightforwardly applied not only to the 
logical qubit space but also the physical qubit space or 
the larger space in which the physical qubit space is em- 
bedded in. 



V. SUMMARY AND DISCUSSION 

In summary we have developed a continuous dynami- 
cal decoupling framework from a geometric perspective. 
Within this framework, the need to perform unphysi- 
cal, arbitrarily strong and fast control is eliminated, and 
bang-bang control can be viewed as a unphysical limit of 
the continuous decoupling in which the control Hamilto- 
nian is unbounded. The decoupling condition is equal to 
require that all time-dependent error vectors s 7 average 
to zero at the end of the decoupling cycle. The trajecto- 
ries of the error vectors are steered by the control Hamil- 
tonian. By expressing control Hamiltonian in terms of 
the generators of SU (n) , the geometric picture provides 
an intuitive way to design the decoupling pulse sequences, 
provided that the structure function of the corresponding 
SU(n) is known. Several examples are given to explicitly 
illustrate how to design the decoupling sequences start- 
ing from a given system-bath Hamiltonian. We also show 
that, in stead of the order of the decoupling group, the 
minimal time needed to finish one continuous decoupling 
cycle. The framework is not restricted to physical qubit 
space or logical qubit space. The ancilla levels or other 
relevant nearby levels can be naturally included in the 
analysis. This is important especially when ancilla levels 
are needed to perform quantum operation or when the 
nearby levels can not be neglected. 

The paper has focused on designing the decoupling 
pulses which are ideal for quantum memory. Less is ad- 
dressed about how to perform quantum operation and 
decoupling at the same time. The continuous decou- 
pling framework developed in this work can be extended 
to systematically treat this problem. [l4|. The detail of 
the analysis is beyond the scope of this work and will 
be presented elsewhere. It is also intriguing to discuss 
the issue of optimal control. In a typical bang-bang de- 
coupling framework a decoupling group is identified, as- 
suming that the corresponding time evolution operators 
can be implemented. However in real system the possi- 
ble control Hamiltonian at our disposal might be limited 
and the desired time evolution operators might not be 
achievable. It is thus nature to ask the complementary 
questions: Giving a set of possible control Hamiltonian, 
what is the optimal decoupling sequence? The geomet- 
ric continuous decoupling framework developed here is 
suitable to answer this question. As a simplest exam- 
ple consider the case where some of the generators of 
SU(n) can not be used as the control Hamiltonian. In 
this case one can use the remaining generators and their 
corresponding rotations to construct a (sub)-optimal de- 
coupling sequence. 

Recently there has been efforts to unify the dynami- 
cal decoupling and the quantum Zeno effect. 0] Roughly 
speaking bang-bang control and quantum Zero effect 
both require a strong interaction with a quantum system. 
This idea shared by two scheme leads to the possible uni- 
fication. It is thus very interesting to see if the continuous 
decoupling scheme in which the arbitrarily strong pulses 
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are not necessary can still be connected to the quantum 
Zcno effect. 



APPENDIX: SU(n) AND SO(N) 

In this appendix we explicitly calculate the correspon- 
dence between the generators of 5(7 (n) and the natu- 
ral generators of SO(N) where N — n 2 — 1 using the 
procedure outlined in Sec IIIII We adapt the following 
convention for the natural generators of SO(N): 



Using this set of generators one can verify the following 
relations 



\L 



(A.l) 



[Ai,A] 
[A 2 ,A] 
[A 3 ,A] 
[A 4 ,A] 
[As, A] 
[A 6 ,A] 
[A 7 ,A] 
[Ag, A] 



where fi < v. Note that the choise of the generators of 
SU(n) is not unique. Unitary transformation on a set of 
generators results in another set of generators. However 
this only amounts to an rotation on error vectors and the 
basis vectors of SO(N). 



1. 517(2) and SO (3) 

Let Ai — a x , A2 = a y , and A3 = a z be the 3 generators 
of SU(2). They satisfy the trace-orthogonality 



+4X 23 + 2X 47 + 2X 65 )A, (A.10) 
-4Li 3 + 2X 46 + 2X 57 )A, 
+AL21 + 2X 45 - 2X 76 )A, 
—2^62 — 2L71 — 2L53 — 2\/3L 58 )\, 
+2Li 6 + 2L 34 — 2X 72 — 2v / 3£s4)A, 
-2X 24 + 2X 37 + 2X 51 + 2V3i 78 )A, 
+2X i4 + 2X 25 - 2X 63 - 2VSL 86 )X, 
+2V3L 45 + 2V3X 67 )A. 

3. 5(7(4) and SO (15) 



Tr(AjAj) = 2Sij. 
It is easy to verify that 



(A.2) 



For SU(n) where n is a power of 2, it is convenient to 
use the product of Pauli matrices to form the generator 
of SU(n). We hence use the following assignment for the 
generators of 517(4): 

Ai = ax ® 1, A 2 = <r 2 <S> 1, A 3 = o- 3 ® 1, (A.ll) 
A4 = 1 ® ax , A5 = 1 (g) (T2, A6 = 1 ® o- 3 , 
A 7 = d ® ax , A 8 = ax <8> <?i , A 9 = ax <8> cr 3 , 
A10 = 0"2 ® Ci, An = ^2 ® ^2, A12 = cr 2 ® CT 3 , 
A13 = cr 3 (g) (Ji, A14 = cr 3 ® cr 2 , A15 = 03 ® cr 3 . 

(A.12) 



A 1; A] = 


+4L23A, 


(A.3) 




A 2 ,A] = 


-4ii 3 A, 


(A.4) 


The generators satisfy the trace-orthogonality 


A3, A] = 


+4Li 2 A. 


(A.5) 


Tr(\ i \ j )=46 ij . 



2. 5(7(3) and SO (8) 

A standard set of generators of the 5(7(3) are the Gell- 
Mann matrices: 



(A.13) 

By matrix manipulation one can verify that the follow- 
ing relations holds 



Ai 




-i 



1 



, A 2 = i ,A 3 = -1 








(A.6) 



A 4 = 



1 

) A 5 = 

1 




1 Ar — 



A 7 




1 



1 
,A 8 = -= I 1 
v 3 \ -2 




They satisfy the trace-orthogonality 
Tr(AiAj) = 2Sij. 



[A X ,A] = 
[A 2 ,A] = 
[A 3 ,A] = 
[A 4 ,A] = 
[A 5 ,A] = 
[A 6 ,A] = 
[A 7 ,A] = 
[A 8 , A], 
[A 9 ,A] = 
[Aio,A] 
[An, A] 
[Ai2,A] 
[Ai3,A] 
[An,A] 
[Ais.A] 



+8^2,3 + 8Lio,i3 + 8Lxx,x4 + 8L1245) A 
—81/1,3 — 8L 7j i3 — 8L8,i4 — 8X9,15) A 
+8X1,2 + 8L 7 ,io + 8Lg,ii + 8X9,12) A 

+8X5,6 + 8Xg,g + 8X11,12 + 8X14,15) A 

—8X4,6 — 8X 7 ,g — 8X10,12 — 8X13,15) A 
+8X4,5 + 8X 7 ,g + 8X10,11 + 8X13,14) A 
+8X2,13 — 8X 3 ,io + 8X5,9 — 8X6,8) A 
—8X2,14 — 8X 3 ,n — 8X4,9 + 8X6, 7) A 
+8X2,15 — 8X3,12 + 8X4,8 — 8X5,7) A 
—8X1,13 + 8X3,7 + 8X5,12 — 8X6,n) A 
—8X1,14 + 8X3,8 — 8X4,12 + 8X6,10) A 
—8X1,15 + 8X3,9 + 8X4,11 — 8X5,10) A 
+8X1,10 — 8X2,7 + 8X5,15 — 8X6,14) A 
+8Xi,n — 8X2,8 — 8X4,15 + 8X6,13) A 
+8X1,12 — 8X2,9 + 8X4,14 — 8X5,13) A 



(A.9) 
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